Single-shot initialization of electron spin in a quantum dot using a short optical pulse 
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We propose a technique to initialize an electron spin in a semiconductor quantum dot with a 
single short optical pulse. It relies on the fast depletion of the initial spin state followed by a pref- 
erential, Purcell-accelerated desexcitation towards the desired state thanks to a micropillar cavity. 
We theoretically discuss the limits on initialization rate and fidelity, and derive the pulse area for 
optimal initialization. We show that spin initialization is possible using a single optical pulse down 
to a few tens of picoseconds wide. 



Semiconductor quantum dots (QDs) constitute good 
candidates for the implementation of quantum informa- 
tion processing in the solid state. In the last ten years, 
the coupling of a single QD transition to an optical cav- 
ity mode has allowed numerous cavity quantum electro- 
dynamics realizations, both in the weak and strong cou- 
pling regimes pHa ■ To benefit from the microsecond spin 
coherence time [6[ , for an electron trapped in a quantum 
dot, spin qubit i nitia lization 

0-0, manipulation [Tol [TT|| 
and readout [H, [l3| are also being widely explored. At 
the interface, several proposals have focused on the use of 
light-matter coupling to achieve remote qubit entangle- 
15j . Here the potential of cavity-QD systems 
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for short-pulse spin initialization is explored. 

High-fidelity electron spin state preparation has been 
obtained by several groups worldwide, using different 
configurations: magnetic field oriented parallel or per- 
pendicular 0, ||| to the optical beam for electron spin 
initialization. Only the perpendicular (Voigt) geometry 
allows initialization times shorter that the typical spin 
coherence time: as shown by Emary et al an initial- 
ization fidelity of 99% can be reached in approximately 
5 nanoseconds in such a configuration. In this Letter we 
show that spin initialization can be performed with a sin- 
gle optical pulse down to a few tens of picoseconds wide, 
using a QD coupled to a high quality factor (Q) optical 
cavity mode, such as the QD-micropillar system sketched 
in Fig. la. 

The corresponding four-level system, in Voigt geom- 
etry, is shown in Figure lb: the low-energy levels are 
the Zeeman-split electron states with spins along the x 
axis, that is |e ± ) = (|t) ± ||)) /V%, where |t) and ||) de- 
note the electron spin states quantized along the z axis. 
The high-energy levels are the Zeeman-split trion states 
I**) = (It) II) - II) It»(l1r> ± W))A where and |4) 
denote the heavy-hole spin states in the z direction. In 
the case of pure heavy-hole states, the \t + ) —> |e + ) and 
\t~) — > \e~) transitions are X-polarized (linearly polar- 
ized parallel to the magnetic field), while the \t + ) — > |e _ ) 
and \t~) — > |e + ) transitions are Y-polarized (perpendicu- 
lar to the magnetic field) [§| . We consider a cavity mode 
in resonance with the |i + ) — > \e~) transition, accelerating 



a) 



laser 

* 



< 



C) 3 



cn 2 



b) 



\t + ) 



a 



r<F p+ D x / 



o 







I 


F P=° 








\ F P =30 




y 









1.5 



1.0 



0.5 



\e + ) 



|e"> 



5 10 15 20 25 



Q 



Vr 



FIG. 1. a) QD-micropillar device in Voigt geometry, b) 
Four-level system and corresponding optical transitions (solid 
line: X-polarized, dashed line: Y-polarized), with Q and F the 
pumping and relaxation rates, respectively, c) Initialization 
time To for F p =0, 10 and 30 as a function of fi/r. 



emission in the mode by a Purcell factor F p . 

The initialization process described in Rcf. Q relies 
on the resonant excitation of the \t + ) —> \e + ) transition 
by an X-polarizcd laser: the system can then relax either 
in the desired state |e _ ), or return in |e + ) where it can 
be excited once again. The population remaining in |e + ) 
thus decreases exponentially with a characteristic time r, 
given in this case by p where T is the spontaneous emis- 
sion rate of the electron-trion transitions. Near-unity fi- 
delity initialization can thus be obtained using optical 
pulses of durations At such that At >• r. 

Here we show that the selective Purcell acceleration 
of the \t + ) — > |e~) transition, not only decreases the 
value of this characteristic time down to r = ,„ ttyp , 
but also allows to initialize with optical pulses of dura- 
tions At shorter than r. In this regime the initialization 
process can be schematically viewed as a two-step pro- 
cedure: complete depletion of |e + ) towards |t + ) using a 
pulse of carefully controlled area, followed by a preferen- 
tial, Purcell-accelerated relaxation towards |e _ ). In the 
following we analyze the initialization rate and fidelity, 
and derive the pulse area required for optimal spin ini- 
tialization. We finally show that micropillar- QD systems 
with the required properties can be realized using state- 
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of-the-art fabrication techniques. 

We consider an X-polarized laser field applied in res- 
onance with |e + ) — > \t + ); the second X-polarized tran- 
sition is detuned by an amount T,b from the laser field 
frequency. The laser being non-resonant with the optical 
mode, the coupling with these transitions occurs mainly 
through the pillar leaky modes, so that we can assume an 
equal pumping rate f2 for both of them. We thus use the 
hamiltonian given by Emary et al. except that we 
define Q as the Rabi frequency of the resonantly driven 
transition, which is twice higher than their definition. 

We determine the properties of our system by solving 
the master equation for the density matrix p in the Lin- 
blad form (h = 1 units) [§]. We assume for simplicity 
that the Purcell effect only alters the spectrally-matched 
\t + ) — > |e _ ) transition whose total spontaneous emission 
rate becomes (F p + 1)T (that is, F p T in the cavity mode, 
and r in the micropillar leaky modes), the other decay 
rates remaining equal to T. We thus use the Lindblad 
operators of Rcf. [8[, with the same value V = 1.2 peV 
(r _1 = 0.55ns), except for a term (F p + 1)T instead of T 
for the accelerated transition. 

We first consider a continous-wave (CW) pumping rate 
f2 instantaneously switched on at time t = 0. In the 
high-Es limit, the varying terms in the hamiltonian os- 
cillate so rapidly that their contribution is averaged to 
zero, and the density matrix elements Pij{t) tend to- 
wards constant values pf) which are all equal to zero 
1. In the long-time limit, the dif- 



except for p^} e 

ferences Pij(t) — pf^ decrease proportionnaly to e T o, 
where the initialization time To can be determined an- 
alytically. Adapting the calculation of Emary et al [H 
to take into account the Purcell acceleration, we find a 
more general form for Tq as a function of the ratio r 
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The evolution of To is shown in Fig. lc: as expected, 
the initialization time is greatly shortened by the Pur- 
cell acceleration of the \t + ) — > \e~) transition. For high 
enough values of ft (r ^ F p + 1), an asymptotic regime is 
reached where To tends towards the characteristic time 
r = ir xtvt i m this regime the initialization speed is not 
limited by the pumping rate but solely by the desexci- 
tation rate of \t + ) towards |e — ). The time required to 
obtain 99% initialization fidelity is approximately 5To, 
which is at best 10r _1 in the case F p = (no cavity), 
but decreases to 0.25r _1 when F p = 40. 

If we now reckon with the complete time-dependant 
hamiltonien (finite Ss), we find that an undesired pop- 
ulation will always remain in the state |e + ). This arises 
from the residual Rabi oscillation in the non-resonantly 
driven, X-polarized transition from |e~) to \t~). This os- 
cillation, with a generalized Rabi frequency y/TP 
and an amplitude proportionnal to -==== , allows a re 
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FIG. 2. a) Population of e J) as a function of time, numer- 
ically integrated for F p = 30 and B — 6 T, using different 

pumping rates around the optimal value Q. op t = 2 • '-') 
Final initialization error as a function of the magnetic field, 
for F p — 30, 40 and 50, for the optimal pumping rate. 



laxation from \t~) to |e + ). To determine the correspond- 
ing populations, we remark that in the long-time limit 
the density matrix elements pij (t) still tend towards con- 
stants pfj\ except for p e - t - which tends towards ge~ l ^ bt , 
with g a complex number. Injecting these matrix ele- 
ments into the master equation leads to the determina- 



tion of the asymptotic values for the populations p. 
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with the denominator D given by: 



ID (1) 
(F p + 3)n 2 ] ID (2) 
(3) 
(4) 



D = (8 + 8Fp + Fp 2 )r 2 + 4(F p + l)Z 2 B + 2(Fp + 2)n 2 (5) 

The residual population in |e + ) is an increasing func- 
tion of the pumping rate f2; thus, once the latter is set 
to an appropriate value for fast spin initialization, ris- 
ing it becomes counterproductive. We find that choos- 
ing Q op t = - = (- F p+ 1 ) r as the optimal pumping rate 
provides a good compromise between initialization speed 
and fidelity. Fig. 2a shows the evolution of p e - e -(t) for 
various pumping rates, with F p = 30 and B = 6T mag- 
netic field, with Yjb determined using the electron and 
hole transverse g factors g\ = —0.46 and g\ = —0.29 0. 
A pumping rate significantly lower than tt op t slows the 
initialization rate; a higher pumping rate enhances the 
residual Rabi oscillations in the \e~) — > \t~) transition. 

We now calculate the corresponding initialization er- 
ror, assuming the CW laser field is turned off at time 
t = At with At >> t, so that at the end of the pulse the 
populations have already converged towards their asymp- 
totic values p^ . At the end of the optical pulse, the 
population remaining in the state relaxes towards 



|e + ) with probability 
|e~) with probability 



-Fp+2 
-Fp + 1 
F p +2 



and towards the desired state 
Similarly, the population in 
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the trion state \t~) relaxes towards |e + ) and |e~) with 
probabilities |. Thus the initialization error we consider 
is the residual population in |e + ) after the pulse, given 
bv o (oo) - o (0) H l —o 



(0) 



We denote e 



this final initialization error, which is found to be: 



2(F P + 2) 2 r 2 
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In Fig. 2b, e is plotted as a function of the magnetic field 
for various Purcell factors, using the optimal pumping 
rate ft op t- The error increases when F p increases, due to 
the higher value of £l pt- As expected, the obvious way 
to limit the initialization error is to increase £& in the 
denominator D, and thus the magnetic field. 

A major feature of the present initialization procedure, 
however, is that the practical implementation of spin ini- 
tialization does not require a pulse duration At higher 
than t. A short optical pulse of carefully controlled area 
can indeed transfer the population from the |e + ) state to 
the \t + ) state; in a second step, the preferential, Purcell- 
accelerated relaxation of \t + ) leads to a nearly complete 
initialization in |e~). For this method to work, a strong 
Purcell factor is needed, so that the transition back to 
|e + ) is quenched as much as possible. The second issue 
is that, for the unintended |e~) — > \t~) transition to re- 
main unaffected by the laser pulse, must be much 
higher than the laser spectral width. 

To calculate the required pulse area for optimal ini- 
tialization, we consider the case of a square optical pulse 
with a constant pumping rate H, in the limit F p >> 1 
and Eg much greater than both and -. In such a case, 
with the laser field switched on at t = 0, the population 
of the |e + ) state simply evolves as : 



Pe+e+ (t) 
Pe+e+ (0) 
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D, 2 \ sin At 
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(7) 



with A = y fl 2 — (i) . We remark that p e + e + reaches 
zero at the end of the optical pulse if the latter's duration 
and pumping rate are related by: 



arctan At 



XAt 



(8) 



The corresponding pulse area, 9 = ft At, can be deduced 
by numerical solving of Eq. (|8|). 

Fig. 3a shows this pulse area for various Purcell fac- 
tors, as a function of the pulse duration. The limit 
At << t leads to a pulse area 9 = it for complete deple- 
tion: this is a well-known property of two-level systems, 
valid as long as the pulse duration is negligible with re- 
spect to the relaxation times. For longer pulses, however, 
the required area increases with At. When At >> t an 
asymptotic behaviour is obtained where 9 ps — . This 
corresponds to O ~ •=■, which is the above- mentionned 
optimal pumping rate fi D pt 
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FIG. 3. a) Pulse area for optimal initialization for F p = 30, 
40, and 50, as a function of the square pulse duration, as 
deduced from Eq. |[HJ. b) Initialization error corresponding 
to this pulse area, for B = 6T, B = 10T, and in the limit of 
infinite magnetic field, with a Purcell Factor F p — 40. 



If this optimal pulse area is used, the limit F p » 1 en- 
sures that the population excited in \t + ) relaxes towards 
|e + ) only, leading to a perfect spin initialization. In prac- 
tice, however, the finite values of the Purcell factor and 
magnetic field lead to a non-zero initialization error, that 
we investigate numerically. Fig. 3b plots the initializa- 
tion error as a function of pulse duration, for 6T, 10T, 
and infinite magnetic fields, using F p = 40. As expected, 
for pulse durations At » t the errors converge towards 
the values described in Fig. 2b, as given by Eq. ((6|) for 
q _ i _ (F p +i)r ^ rp ne ogcjrja^oj-y behaviour arises from 
the non-adiabatic switching of the laser field at t = and 
t = At: as shown in Fig. 2a, this leads to residual Rabi 
oscillations in the |e _ ) — > \t~) transition. In the limit of 
a very short pulse and for an infinite magnetic field, the 
initialization error tends towards: 
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F p + 2' 



(9) 



(F P +i)r 



which is the population transferred from |e + ) to \t + ) mul- 
tiplied by the probability of the unwanted relaxation from 
\t + ) back to |e + ). However, in this regime, and for realis- 
tic magnetic fields, the condition E^Ai >> 1 is no longer 
satisfied: the high spectral width of the pumping laser 
precludes any efficient initialization, as the |e~) \t~) 
transition is strongly excited. 

For these short optical pulses, the approximation of a 
square pulse is often inappropriate; yet the validity of the 
method does not depend on the pulse shape as long as 
the pulse area 8 = J Q.(t)dt is carefully controlled. Fig. 
4a shows the numerical calculation of the |e + ) and |e~) 
populations for a 6 T magnetic field and a F p = 40 Pur- 
cell factor, using a gaussian pulse with a 40 ps (0.07r _1 ) 
full-width at half-maximum (FWHM) . In such a case the 
magnetic field is high enough to ensure that the unin- 
tended |e~) to transition is non-resonant with the 
laser pulse, which results in the absence of depletion of 
the |e~) state during the pulse. The preferential relax- 
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FIG. 4. a) Optical pulse normalized intensity and correspond- 
ing populations of |e+) and |e— ) states, as a function of time 
(40 ps FWHM gaussian pulse, with F p = 40 and B = 6T). b) 
Initialization error as a function of the gaussian pulse width 
for 6 T, 10 T, and infinite magnetic fields, with F p — 40. 
Insert: required pulse area for optimal initialization using a 
gaussian pulse, for Fp — 30, 40, and 50. 



ation towards |e~) leads to a final initialization fidelity 
around 99%, measured at time t = 140ps (0.25r _1 ); 
the initialization process is thus approximately 40 times 
faster than without a cavity. To achieve such a fidelity 
with a 20 ps FWHM pulse, a Purcell factor F p = 50 and 
a magnetic field B = 10 T are required. 

Fig. 4b shows the expected initialization error for 6T, 
10T, and infinite magnetic fields, using F p = 40. As ex- 
pected, the initialization error in the limit of infinite mag- 
netic field and very short pulse duration remains given 
by Eq. [9j The main difference with Fig. 3b is that the 
smooth evolution of fl(t) allows a monotonic decrease of 
the initialization error when the pulse width increases. 
As long as the pulse FWHM is larger than a few tens of 
picoseconds, thanks to the limited spectral width of the 
gaussian pulse, the finite value of the magnetic field does 
not significantly alter the initialization efficiency. 

The insert of Fig. 4b summarizes the required pulse 
area, as a function of the pulse FWHM, displaying the 
same behavior than in Fig. 3a for square pulses: the 
pulse area for optimal initialization grows above it when 
the pulse width increases. We note that the pulse area 
shows no significant dependence on the magnetic field in 
the 6T— 10T range, and that no fine-tuning is required for 
practical implementation: for example, a 2% shift of the 
pulse area will result in an increase of the initialization 
error of the order of a few 10 -3 . 

We now show that this spin initialization procedure 
can be implemented in realistic QD-cavity systems. The 
current state-of-the-art of nano-fabrication techniques al- 
lows to realize, in a both deterministic and scalable 
way, QD-micropillar systems whose cavity mode is spa- 
tially and spectrally matched with a selected QD transi- 
tion [lf|. Micron-size pillars can be realized with qual- 
ity factors up to 10 5 



V w 1 fim 3 and effective index n w 3.4) with a qual- 
ity factor Q = 2.5 10 , embedding a single InAs/GaAs 
QD with / = 10 oscillator strength (transition wavelenth 
A = 0.94/im) will present a Purcell Factor F p w 40 while 
staying in the weak-coupling regime 0]. 

A key condition to the validity of this initialization pro- 
cedure is that the undesired relaxation from \t + ) back to 
|e + ) does not experience any significant Purcell acceler- 
ation. Using the transverse g factors of Ref. [|[ , we find 
that a 6T magnetic field corresponds to a 160 fieV de- 
tuning between the \t + ) — > \e + ) transition and the cavity 
mode, whereas the mode resonance half-width is 25/ieV 
for Q = 2.5 10 4 : this corresponds to a Purcell factor 
F p rj 1 (F^ 0.3 for B = 10 T). The small Purcell 
effect experienced by the Y-polarized \t~) — > |e + ) transi- 
tion (the Zeeman splitting between the Y-polarized tran- 
sitions is only 60/^eU for B = 6 T) does not alter the 
validity of the method. 

We also point out that the presence of heavy /light-hole 
mixing in the QD can affect the trion states, resultingin 
a rotation of the polarization axis in the x — y plane [9j ; 
furthermore, a small alteration of the beam polarization 
may also occur in the Bragg mirrors, due to the resid- 
ual ellipticity of the micropillar (cllipticities of the order 
of 5.10~ 4 are the current state-of-the-art For sim- 

plicity the orthogonal polarization states of the optical 
transitions have been called X and Y , whether or not 
they correspond to linear polarizations along the x and 
y axes of Fig. la. For efficient spin initialization the 
laser polarization has to be adjusted to ensure negligible 
excitation of the Y-polarized transitions. 

Having derived the optical pulse area required for op- 
timal initialization, as well as the corresponding fidelity 
and experimental conditions, we conclude that spin ini- 
tialization can be implemented using short optical pulses 
down to a few tens of picoseconds wide. This work is par- 
tially supported by ANR project n°ANR-JCJC-09-0136. 



17] . For example, a 3fim di- 



ameter AlAs/GaAs micropillar (effective mode volume 
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